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Abstract 

The  research  project  studied  the  benefits  of  Magnetogasdynamic  (MGD)  ac¬ 
celeration  on  a  scramjet  nozzle.  MGD  is  a  technology  that  relies  on  electromagnetic 
fields  to  extract  and/or  add  energy  to  flow  fields.  The  scramjet  engine  in  the  “AJAX” 
concept  proposed  by  D.I.Brichkii  et  al.  (2001)  of  St.  Petersburg,  Russia  would  utilize 
an  MGD  power  generator  in  the  diffuser  which  could  potentially  provide  electrical 
power  for  the  aircraft  support  systems  and  an  MGD  accelerator  in  the  nozzle  of  the 
engine.  Using  an  AFRL/VAAC  CFD  code  that  was  modified  for  MGD  computations, 
this  project  investigated  the  effect  conductivity  and  load  factor  had  on  the  specific 
thrust  and  efficiency  of  an  MGD  accelerator  composed  of  segmented  electrodes  in 
the  nozzle  of  a  scramjet  engine.  For  a  load  factor  of  1.6,  results  showed  a  95%  in¬ 
crease  in  specific  thrust.  When  MGD  generation  was  employed  at  a  load  factor  of  0.8, 
a  31%  decrease  in  specific  thrust  was  calculated.  The  solutions  included  high  wall 
temperatures  which  need  more  investigation  for  MGD  acceleration  to  be  practical. 
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STUDY  OF  MAGNETOGASDYNAMIC  FLOW 
ACCELERATION  IN  A  SCRAMJET  NOZZLE 

1.  Introduction 

Since  the  Wright  brother’s  first  self  powered  flight  on  December  17,  1903 [16], 
humans  have  achieved  great  advances  in  aviation.  From  the  supersonic  X-l  flown  by 
Chuck  Yeager  on  October  14,  1947[16]  to  Burt  Rutan’s  2005  X-Prize  winning  Space 
Ship  One,  the  world  is  quickly  getting  smaller.  The  transportation  industry  counts 
down  the  days  until  a  flight  from  New  York  to  Tokyo  only  takes  two  hours [7,  2],  In 
this  industry,  speed  is  addictive. 

Rockets  and  Ramjets  have  set  the  pace  at  Mach  6.  However,  these  propulsion 
systems  are  limited  by  weight  and  drag;  and  new  technologies  are  being  investigated. 
In  the  designs  prior  to  the  X-l,  the  greatest  limitation  was  structural.  With  airplanes 
falling  apart  as  pilots  attempted  to  break  through  the  sound  barrier,  the  process 
was  primarily  trial  and  error[16,  33].  At  the  expense  of  labor  and  lives,  earlier  ad¬ 
vancements  required  deep  pockets  and  long  turnaround  times.  Today’s  engineering 
methods  include  ground  testing  and  computational  fluid  dynamics  (CFD)  to  explore 
new  frontiers. 

Sustained  hypersonic  flight  has  been  that  new  frontier  for  over  fifty  years[16]. 
Current  hypersonic  propulsion  systems  include  rocket  and  ramjet  engines.  However 
these  systems  have  exhausted  their  speed  potential.  Rockets,  though  less  susceptible 
to  drag,  are  held  back  by  the  weight  of  the  oxidizers  that  they  must  carry.  Airbreathing 
ramjet  propulsion  eliminates  the  need  for  internal  oxidizers. 

Ramjets  operate  by  slowing  down  the  hypersonic  flow  to  subsonic  mach  num¬ 
bers  for  combustion [17].  Atmospheric  oxygen  is  mixed  with  fuel  and  ignited.  How¬ 
ever,  engine  unstart  can  occur  when  too  much  fuel  is  introduced  to  the  mixture. 
This  increases  the  back  pressure  expelling  the  shock  train  from  the  inlet.  Even  with 
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proper  fuel  flow  control,  ramjets  are  not  always  effective.  When  approaching  Mach 
6,  the  process  of  slowing  down  the  flow  field  creates  too  much  drag  for  the  thrust  to 
overcome[24,  16,  2],  This  is  essentially  the  operational  limitation  that  leads  to  the 
development  of  scramjets. 

1.1  Scramjet  Research 

Supersonic  combustion  ramjets  (scramjets)  eliminate  the  limitations  that  are 
inherent  in  rockets  and  ramjets.  Scramjets  do  not  rely  on  internal  oxidizers,  nor  do 
they  require  subsonic  combustion [17].  However,  there  are  significant  challenges  that 
must  be  met  for  the  sustained  hypersonic  flight  in  the  scramjet.  Examples  include 
engine  unstart,  and  the  fact  that  the  scramjet,  like  the  ramjet,  needs  a  boost  from  an 
external  supersonic  propulsion  system  to  get  started[16,  2], 


Figure  1.1  University  of  Queensland’s  HyShot[19]. 


The  University  of  Queensland’s  HyShot  Flight  Program  is  one  of  several  that 
has  investigated  scramjet  propulsion.  Their  two  flight  dates  were  30  Oct  2001  and 
30  July  2002.  Their  next  launch  is  scheduled  for  November  2005 [19].  A  picture  of 
their  scramjet  engine  and  rocket  stage  is  shown  in  Figure  1.1.  The  rocket  was  used  to 
bring  the  scramjet  up  to  330  km.  The  rocket  and  scramjet  engine  were  then  directed 
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back  down  and  the  second  rocket  stage  was  fired.  The  scramjet  then  separated  and 
propelled  itself  for  5  seconds  at  Mach  7.6  between  35  km  and  23  km  on  hydrogen 
fuel  [19,  16]. 

The  Hyshot  team  had  some  difficulties  unique  to  their  design.  Due  to  their 
vertical  trajectory  during  the  launch,  accelerations  reached  a  peak  of  60g  before  de¬ 
caying  to  30g  which  could  cause  structural  damage  during  the  first  second  of  flight 
[16,  19].  Recent  ground  tests  and  CFD  studies  at  the  University  of  Queensland  show  a 
promising  revival  of  their  design.  Rated  for  Mach  10  conditions  at  30  km  in  altitude, 
the  HyShot  team  promises  to  beat  the  record  of  its  competitor  in  November  2005  [20]. 

NASA  is  The  University  of  Queensland’s  main  competitor.  NASA’s  $250  million 
X-43A  Hyper-X  has  had  two  successful  flights.  Unlike  the  Australian  HyShot,  it 
incorporated  the  engine  into  the  aircraft.  This  is  done  so  that  the  underside  of  the 
fuselage  provides  external  combustion  before  the  inlet  and  external  thrust  surface  after 
the  nozzle.  As  shown  in  Figure  1.2,  the  X-43A  is  propelled  by  a  Pegasus  missile  that 
is  launched  by  a  B-52  bomber.  The  first  flight  occurred  on  March  27,  2004  during 
which  it  reached  Mach  7  at  100,000  ft [16,  34],  On  its  last  flight  it  reached  nearly 
Mach  10,  or  7000  mph,  as  it  flew  at  an  altitude  of  110,000  ft [25].  The  current  design 
has  a  few  shortfalls.  The  X-43A  is  aerodynamically  unstable,  and  the  engine  has  a 
tendency  to  overheat  [16].  This  is  coupled  with  the  control  problem  of  too  much  fuel 
getting  into  the  combustor  which  increases  back  pressure  and  can  eventually  cause 
unstart  [16,  2], 

1 . 2  Magnetogasdynamics 

In  response  to  flow  problems,  some  researchers  have  started  investigating  mag¬ 
netogasdynamics  (MGD)  to  control  scramjet  flow  fields[16,  24,  7].  MGD  is  the  study 
of  electrically  conducting  gases.  At  velocities  below  Mach  9  at  which  sufficient  spon¬ 
taneous  or  thermal  ionization  does  not  occur  [31],  electron  beams  can  be  utilized  to 
ionize  the  flow.  MGD  has  been  investigated  for  a  variety  of  aerospace  applications 
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Dryden  Flight  Research  Center  ED97  43968-04 
B- 52  CARRY-ALL:  This  artist's  concept  depicts  the  Hyper-X 
research  vehicle  riding  on  a  booster  rocket  prior  to  being  launched 
by  Dryden  Flight  Research  Center's  B-52  at  about  40,000  feet. 


Figure  1.2  NASA’s  X-43A[25], 


from  boundary  layer  control  around  around  a  reentry  vehicle  [29]  to  flow  field  acceler¬ 
ation  in  a  nozzle. 

Ohms  law  plays  an  intricate  role  in  accelerating  flow.  The  details  of  the  law 
are  discussed  in  Chapter  2.  As  an  overview,  the  flow  field  can  be  accelerated  with 
the  addition  of  electrical  energy  to  the  flow,  or  it  can  be  decelerated  with  the  gener¬ 
ation  of  electrical  energy  from  the  flow.  One  concept  known  as  “AJAX”  originated 
in  St.  Petersburg  Russia  with  A.  Kuranov  and  E.  Sheikin  at  the  Hypersonics  Sys¬ 
tems  Research  Institute[16]  and  involves  both  flow  held  acceleration  and  generation. 
Figure  1.3  shows  a  diagram  of  the  “AJAX”  concept.  The  Magneto-Plasma-Chemical 
Engine  (MPCE),  also  referred  to  as  a  scramjet  with  MGD  bypass,  was  developed 
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within  its  framework.  The  engine  design  incorporates  an  MGD  generator  that  serves 
two  purposes.  It  provides  energy  for  the  flight  systems  and  an  MGD  accelerator  which 
utilizes  the  Lorentz  force  vector  in  a  direction  that  accelerates  the  flow  field,  and  the 
generator  serves  as  a  means  of  reducing  the  flow  velocity  for  combustion  without 
significant  pressure  drop [16,  7,  14]. 


Figure  1.3  The  “AJAX”  concept[9]. 

The  implementation  of  MGD  is  based  on  several  conditions.  These  include 
method  of  ionization,  and  magnet  and  electrode  placements.  The  flow  can  be  ionized 
by  injection  of  ionized  particles  or  through  use  of  an  electron  beam  ionizer [11].  The 
electron  beam  ionizer  is  typically  suggested  since  it  produces  the  best  current  density 
to  weight  ratio  to  date  [24], 

The  Lorentz  force,  or  j  x  B,  is  the  driving  force  for  MGD  generation  and  accel¬ 
eration  in  the  scramjet.  In  order  for  this  force  to  accelerate  flow  in  the  axial  direction, 
the  magnetic  flux  density,  B ,  and  the  current  density,  j,  must  be  perpendicular  to 
the  direction  of  flow,  u  as  shown  in  Figure  1.4  below. 

1.3  Current  Research 

This  thesis  is  part  of  a  progressive  effort  to  increase  understanding  of  MGD  flow 
control  in  a  scramjet  engine.  The  ultimate  goal  is  to  determine  conditions  required  for 
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a  self-sustained  scramjet  engine  that  utilizes  MGD  generation  to  power  its  onboard 
electronics  with  enough  power  left  over  to  accelerate  the  flow  in  the  nozzle  for  added 
thrust.  This  particular  study  is  focused  on  determining  the  conditions  that  improve 
efficiency. 

Figure  1.5  shows  the  whole  scramjet  engine  as  proposed  by  Gaitonde[14],  This 
study  investigates  an  accelerator  with  three  electrodes  as  shown  in  Figure  1.4.  The 
electrodes  are  the  vertical  plates  on  either  side  of  the  nozzle. 

Harrington’s  work  focused  on  flow  held  interaction  and  what  parameters  caused 
the  greatest  effects  for  continuous  electrodes[16].  His  work  agreed  with  that  of 
Gaitonde  which  stated  that  significant  Joule  heating  effects  at  the  electrodes  made 
efficient  how  held  acceleration  impractical  for  real  implementations[16,  10].  However, 
their  work  also  emphasized  that  an  a  study  of  efficiency  should  be  performed. 

The  objective  of  this  computational  study  was  to  address  efficiency  while  also 
determining  the  effectiveness  of  differing  conductivity  pattern  and  load  factor.  Both 
acceleration  and  generation  cases  were  considered. 

With  Maxwell’s  equations  and  Navier-Stokes  equations  serving  as  the  corner¬ 
stone,  the  governing  equations  for  MGD  are  summarized  in  Chapter  2.  Then,  Chapter 
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Figure  1.5  Simplified  Mach  8  scramjet  model  design  [10] 

3  contains  a  description  of  the  computation  setup  as  well  as  the  specific  thrust  and 
efficiency  models  that  were  employed.  The  results  follow  in  Chapter  4  which  presents 
overall  specific  thrust,  quasi  one  dimensional  interpretations  of  specific  thrust,  and 
three  dimensional  flow  held  plots.  Chapter  4  ends  with  a  discussion  of  the  efficiency  of 
MGD  acceleration.  Chapter  5  concludes  with  remarks  on  what  was  discovered  from 
the  research  and  what  may  be  considered  for  future  research. 
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2.  Governing  Equations 


2.1  Navier-Stokes  Equations 

The  Navier-Stokes  equations  (NSE)  govern  the  behavior  of  fluid  motion.  They 
represent  conservation  of  mass,  momentum,  and  energy.  When  radiative  heat  transfer 
and  body  forces  are  neglected,  these  laws  are  expressed  in  the  following  differential 
equations  [7]: 


dp  d(pVj) 
dt  dxi 


(2.1) 


S(pVi)  d(pViVj)  _  _dp_  dry 

dt  dxj  dxi  dxj 


d(pet)  d(phtVj)  =  dj-qj  +  , 

dt  dxj  dxj  {  J 

where  i  and  j  refer  to  directions  x  and  y. 

Eqns.  2.2  and  Eqn.  2.3  contain  the  stress  tensor,  Ty,  and  the  heat  flux  vector, 
qv  which  are  defined  below, 


Ti 


dVi  dVj 


2X  dVk 
35ijfifdxk 


(2.4) 


Qj 


dT 


(2.5) 


The  equations  above  incorporate  molecular  viscosity,  pf,  and  thermal  conduc¬ 
tivity,  Kf,  which  are  found  using  Sutherland’s  Law  and  the  local  temperature,  T. 
These  equations  are  listed  below, 


To  +  S,j,f 

T  +  S H 


(2.6) 
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Kf  ~  K0,f 


(2.7) 


/  T  \  *  T0  +  SKf 
V To)  T  +  SKf 

where  p,o j  =  1.716  x  10~5,  T0  =  273  K,  S/lf  =  111  K,  n0j  =  2.41  x  10-2,  and  SKf  =  194 
K  for  air [7]. 

2.1.1  Ideal  Gas  Law.  The  NSE  above  have  five  equations  and  six  unknowns, 
(p,u,v,w,p,et).  In  order  to  complete  the  computations,  the  ideal  gas  law  is  used [7]. 
The  ideal  gas  law  assumes  that  air  is  a  calorically  perfect  gas  so  the  pressure  is  directly 
proportional  to  density,  p,  and  temperature,  T, 


p  =  pIZT  (2.8) 

where  IZ  is  the  specific  gas  constant  for  air. 

Total  enthalpy,  ht,  and  total  energy,  et,  result  from  the  ideal  gas  law.  These 
terms  are  defined  bclow[7], 


ht  = 


h^-VM 


(2.9) 


7  P 

et  -  hi - 

P 


(2.10) 


2.1.2  Vector  Form  of  the  Navier- Stokes  Equations. 
form  of  the  NSE  is  shown  below, 


The  conservative  vector 


<9U  <9E  dF  dG  _ 

dt  dx  dy  dz 

U  is  called  the  vector  of  conservative  variables,  given  by  the  following: 


(2.11) 


U  =  [p  pu  pv  pw  pet\ 


T 


(2.12) 


2-2 


E.  F.  and  G  are  the  total  fluxes  in  the  corresponding  x,  y,  and  z-directions.  These 
are  each  separated  into  an  inviscid  flux  vector  minus  a  viscous  flux  vector  [7]: 


E  =  Ei  —  Ev  = 


G  =  G; 


F  =  Fi-  Fv  = 


pu 

pu 2  +  p 
puv 
puw 
phtu 

pv 

puv 

pv2  +  p 
pvw 
phtv 


7~xx 

TXy 

^~xz 

UTXX  +  VTxy  +  WTXZ  +  qx 

0 


r. 


xy 


T, 


yy 


To 


yz 


^TTXy  T  TTyy  -\-  UJTyZ  T  ([y 


Gv  = 


pw 

0 

puw 

Txz 

pvw 

— 

Tyz 

pw2  +  p 

T zz 

phtw 

UTXZ  +  VTyz  +  WTZZ  +  qz 

(2.13) 


(2.14) 


(2.15) 


For  non-MGD  cases,  the  source  term,  S,  is  equal  to  the  zero  vector, 


S 


0 

0 

0 

0 

0 


(2.16) 


2-3 


2. 2  Special  Relativity 

Since  MGD  utilizes  electromagnetic  fields,  which  essentially  consist  of  light 
waves,  special  relativity  was  considered  in  the  development  of  the  MGD  equations. 
However,  the  special  relativity  terms  were  determined  to  be  negligible  and  were  not 
implemented  in  the  code. 

2.2.1  Reference  Frames.  There  are  two  reference  frames  associated  with 
MGD  equations.  One  refers  to  the  coordinate  system  of  the  observer  who  may  be 
in  the  cockpit  of  a  scramjet  driven  aircraft  or  a  scientist  in  the  lab.  This  frame  is 
called  the  lab  frame  because  its  usually  associated  with  the  inertial  reference  frame 
of  a  lab.  The  other  frame  refers  to  the  coordinate  system  of  the  moving  fluid.  This 
frame  is  called  the  rest  frame  because  it  is  associated  with  the  coordinate  system  that 
is  resting  on  a  particle  that  is  moving  with  the  fluid [18]. 

2.2.2  Lorentz  Transformations.  The  Lorentz  transformations  (LT)  consist 
of  the  following  equations  that  transform  from  the  lab  frame  to  the  rest  frame, 


x'  =  (3{x  —  Vt) 

(2.17) 

Vx 

t  =/3{t - T) 

(2.18) 

where  the  primes  refer  to  the  rest  frame;  V  is  velocity  magnitude  of  flow  and  c  is  the 
speed  of  light.  In  the  above  equations,  (3  accounts  for  coordinate  system  warping  due 
to  special  relativity  effects  and  is  defined  below: 


P  = 


(2.19) 


In  order  to  transform  from  a  vector  in  one  reference  frame  to  a  vector  in  another 
reference  frame  moving  at  a  substantially  fast  velocity  the  Silberstein  operator  a,  also 
known  as  the  stretching  operator [18],  is  used. 
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(2.20) 


-1  )V(V-G) 

ll^ll2 

This  operator  utilizes  the  f3  term  and  the  velocity  of  the  rest  frame  with  respect  to 
the  lab  frame.  The  vector  G  is  arbitrary.  For  this  study,  G  could  be  the  Lorentz  force, 
electric  field,  or  the  current  density. 

2.3  Maxwell  Equations 

The  Maxwell  equations  (ME)  are  the  governing  equations  for  electromagnetics. 
They  are  valid  in  both  the  lab  frame  and  the  reference  frame [7]. 


aG  =  G  + 


(0 


v  ■  D  =  pe 

(2.21) 

V-B  =  0 

(2.22) 

dB 

(2.23) 

Vx£  =  -ar 

-  -  d  D 

(2.24) 

Vxff=i+  ~m 

where  D  is  the  electric  flux  density,  B  is  the  magnetic  flux  density,  E  is  the  electric 
held,  H  is  the  magnetic  field,  j  is  the  conduction  current  density,  and  pe  is  the  space 
charge. 

The  ME  consist  of  four  equations  and  seven  unknowns.  Three  more  relations 
are  used  to  balance  the  system.  One  of  which  is  Ohm’s  Law,  and  the  remaining  two 
are  the  constitutive  relations [18]. 

2.3.1  Ohm's  Law.  Ohm’s  law  is  defined  in  terms  of  anisotropic  media, 
which  means  that  the  conductivity,  cr,  depends  on  the  orientation  of  the  particle.  The 
anisotropic  quality  is  expressed  with  the  directional  subscripts  %  and  j.  The  expression 
is  simple  in  the  rest  frame, 
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•/  771/ 

J*  =  VjiEj 


(2.25) 


However,  Ohm’s  law  does  not  transform  into  the  lab  frame  as  easily  as  the  ME. 
Instead,  Ohm’s  law  depends  on  the  Maxwell- Lorentz  transformations  (MLT).  The 
MLT  are  derived  after  several  substitutions  and  partial  differentiations  of  the  ME 
and  LT,  Ohm’s  law,  and  the  constitutive  relations [18].  Briefly,  Ohm’s  law  becomes, 

ji  =  (~)WAEJ  +  (V  X  B)j]}  +  peVi  (2.26) 

a 

since  V2  <<  c2  for  flow  velocities  around  Mach  8,  the  (^)  term  equals  one.  Ohm’s 
law  then  becomes, 


ji  —  aji[Ej  +  iy  x  B)j]  +  peVi  (2.27) 

2.3.2  Constitutive  Relations.  The  two  constitutive  relations  for  linear 
isotropic  dielectrics  and  magnetic  materials  expressed  in  the  rest  frame  are  listed 
below,  [18]: 


ee 


B'  =  peH' 


(2.28) 

(2.29) 


where  subscript  0  represents  the  freespace  value. 

With  application  of  MLT,  the  above  constitutive  relations  can  be  transformed 
to  the  lab  frame.  When  V 2  «  c2  and  the  fluid  is  isotropic,  the  constitutive  relations 
become  [18], 


D  =  ee{E  +  (1  -  -T-)!?  x  B] 

eo,e  M Q,e 


(2.30) 
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B  =  ne[H-{  1 


1 


(2.31) 


€e  Me 
eO,e  M 0,e 


)y  x  d\ 


2.3.3  Magnetogas  dynamic  Assumptions.  A  current  moving  through  a  sta¬ 
tionary  object  in  an  electric  field  produces  a  Lorentz  force  in  a  direction  perpendicular 
to  the  field  and  the  direction  of  current  flow.  In  fluids,  a  force  acting  on  the  flow  is 
called  a  body  force.  The  Lorentz  body  force  density  for  anisotropic  moving  media  in 
the  rest  frame  is  defined  below[18], 


fe  =  peE  +  jxB-  | E2Wk  -  ^H2Wnm  +  | V(£2^p)  +  ^\7(H2^-p)  (2.32) 

where  the  subscript  m  represents  the  magnetostriction  term. 

The  above  equation  is  complex,  and  is  only  valid  in  the  rest  frame.  How¬ 
ever,  the  body  forces  can  be  simplified,  for  the  lab  frame,  given  the  following  six 
assumptions  [10]: 

MGD  Assumption  1:  \V\2  «  c2.  The  magnitude  of  the  velocities  in  MGD  are 
much  less  than  the  speed  of  light  so  (3  approaches  one. 

MGD  Assumption  2:  E  &  0(V  x  B).  The  applied  magnetic  field  is  much  greater 
than  the  induced  magnetic  field. 

MGD  Assumption  3:  ^  ~  0,  so  ^  ~  0,  and  VxB  =  pej.  The  electric  field  has 
no  high  frequencies. 

MGD  Assumption  4:  eE2  «  — .  Magnetic  energy  is  orders  of  magnitude  greater 

Me 

than  the  electric  energy. 

MGD  Assumption  5:  j  =  a(E  +  V  x  B).  The  conductivity  is  independent  of 
electric  field,  and  constant  with  frequency.  This  means  that  j'  =  j. 
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MGD  Assumption  6:  /  =  peE  +  j  x  B.  The  body  force  density  is  homogenous 
and  independent  of  electrostriction  and  magnetostriction.  Since  peE  «  j  x  B, 
then  /  =  j  x  B  is  used. 

2.3.4  Maxwell’s  Equations  for  Magnetogas  dynamic  Flow.  The  six  assump¬ 
tions  listed  above  can  be  applied  to  the  original  ME.  In  the  lab  frame  they  are  the 
following[18]: 


-*  BE 

Vx£  =  -a 

(2.33) 

VXB  =  pej 

(2.34) 

V-J  =  0 

(2.35) 

V-B  =  0 

(2.36) 

Ohm’s  law  is, 

j  =  a(E  +  VxB) 


(2.37) 


All  subsequent  equations  are  written  in  the  lab  frame. 

2.4  Magnetogasdynamic  Equations 

The  magnetogasdynamic  equations  incorporate  ME  and  the  NSE  for  MGD. 
These  equations  describe  the  interaction  between  electromagnetic  fields  and  electri¬ 
cally  conducting  gases  in  a  continuum  governed  by  the  MGD  assumptions.  The 
magnetogasdynamic  equations  are  the  same  as  those  given  in  Eqns.  2.11  except  the 
term  S  is  now  given  by: 
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jyBz  jzBy 

jzBx  jxBz 


(2.38) 


ixBy  jyB 


X 


Exjx  Byjy  ~t~  Ezjz 


The  three  middle  rows  refer  to  the  three  components  of  the  (j  x  B)  Lorentz  force. 
The  last  row,  equivalent  to  E  ■  j,  is  the  resulting  energy  interaction  which  accounts 
for  Joule  heating  and  the  work  of  the  Lorentz  body  force. 


2-4-1  Non-dimensionalizing  the  Magnetogasdynamic  Equations.  In  order  to 
avoid  machine  limitations  in  the  computational  code,  the  MGD  equations  are  non- 
dimensionalized.  In  the  following  relations,  dimensionless  quantities  are  flagged  with 
an  asterisk  [7]: 


L* 

P* 

E* 


L 

L0 

P_ 

Po 

E 

E~o 


V* 


B* 


V 

V0 

T 

T0 

B 

Bo 


(2.39) 


These  initial  dimensional  freestream  quantities  with  the  subscript  0  are  used  to 
non-dimensionalize  the  other  variables,  and  all  dimensionless  variables  are  substituted 
into  the  equations.  Several  non-dimensional  parameters  are  formed  during  this  process 
which  help  to  characterize  MGD  flows.  The  Reynolds  number  is  the  ratio  of  inertial 
forces  to  viscous  forces, 


where. 


Re 


VqLq 

vf 


(2.40) 


2-9 


Uf  =  ^  (2.41) 

P 

The  magnetic  Reynolds  number  is  a  measure  of  the  magnitude  of  the  induced  magnetic 
field  compared  to  the  total  magnetic  field: 


Rrn  =  V0L0(ToPo,e  (2.42) 

The  Mach  number  and  magnetic  Mach  number  are  listed  below, 


Go 


Mm.= 


V()\/P0P0,e 


Bn 


(2.43) 

(2.44) 


The  Prandtl  number  is  the  ratio  of  kinematic  viscosity  to  thermal  diffusivity[7]: 


Pr  =  (2.45) 

K, 

The  magnetic  Prandtl  number  is  the  ratio  of  viscous  diffusion  to  magnetic  diffusion: 


Pm  =  <J0vfn  o>e  =  (2.46) 

1x6 

The  interaction  parameter,  Q,  is  used  when  Rrn  <<  1.  Q  is  the  ratio  of  the 
ponderomotive  force  to  the  inertial  force,  which  should  be  about  one.  An  increase  in 
Q  causes  an  increase  in  the  strength  of  the  magnetic  field. 


ctqBqLo  Rrn 

PoV0  =Ml 


(2.47) 


The  ratio  of  the  electric  field  to  the  induced  effects  is  known  as  the  load  factor, 
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-Eo 

B0Vo 


(2.48) 


K  = 


where  /C  is  proportional  to  the  applied  electric  held,  E0,  and  voltage  differential,  A 0. 

In  the  non-dimensional  MGD  equations,  the  non-dimensional  Lorentz  force  is 
Q(f  x  B*)  and  the  non-dimensional  energy  interaction  is  Q(E*  ■  j*). 
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3.  Computational  Set-Up 


3.1  MGD  3-D  Computational  Code 

The  computational  study  employed  a  3-D  non-dimensional  CFD  code.  The 
code  was  written  by  AFRL/VAAC[10,  11,  12,  13,  14,  15,  27,  29].  The  code  solves  the 
Euler  and  Navier  Stokes  equations  with  the  source  term  for  MGD  Lorentz  Forces  and 
Energy  Interaction. 

In  the  Hypersonic  case,  high  speed  fluids  involve  discontinuities  in  the  flow  solu¬ 
tion.  In  order  to  capture  these  discontinuities,  flux  vectors  are  calculated  using  Roe’s 
flux  difference  splitting  scheme  with  van  Leer  harmonic  limited  Monotone  upstream- 
centered  schemes  for  conservation  laws  (MUSCL)[11,  12,  13,  14,  28,  29]  which  provide 
up  to  third-order  spatial  accuracy.  Since  the  Roe  scheme  is  capable  of  producing 
non-physical  expansion  shocks  under  certain  constraints,  an  entropy  fix  was  added  to 
negate  this  effect. 

3.2  Thrust  Calculation  Model 

It  is  common  to  evaluate  the  stream  thrust  function  at  the  inlet  and  outlet  of 
a  nozzle  to  determine  the  thrust  provided  by  the  nozzle.  The  equation  for  specific 
thrust  is  defined  in  terms  of  the  stream  thrust  function, 

F 

~T~  Saexit  (3-1) 

m 

However,  a  better  measure  of  the  effect  each  segment  of  electrode  and  insulated 
region  has  on  the  thrust  was  found  by  plotting  out  the  stream  thrust  function  over 
the  entire  length  of  the  nozzle.  This  way  stream  thrust  could  be  interpreted  based  on 
the  location  within  the  stream. 

The  stream  thrust  function  is  defined  [17]  as, 

Sa  =  —  (3.2) 

m 
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where  the  axial  force  function,  /,  is, 

/  =  pA{l  +  7  M2)  (3.3) 

/  measures  the  axial  force  between  two  stations.  The  above  equation  is  based  on  the 
following  assumptions [16]. 

Sa  Assumption  1:  The  flow  is  steady  or  quasi  steady. 

Sa  Assumption  2:  The  ratio  of  momentum  flux  to  pressure  forces  is  large. 

Sa  Assumption  3:  Thermochemistry  of  the  flow  must  be  known. 

Sa  Assumption  4:  The  velocity  component  at  each  stage  aligned  with  the  thrust  or 
axial  direction  and  perpendicular  to  the  through-flow  area  is  used  to  calculate 
the  Mach  number. 

Steady  flow  inertial  terms  are  large  in  hypersonic  flow,  so  Assumption  1  is 
valid.  Assumption  3  assumes  that  thermochemistry  was  known.  This  can  be  assumed 
since  the  atmosphere  is  primarily  composed  of  N2  and  02  molecules  which  do  not 
significantly  ionize  at  Mach  numbers  less  than  9  [31]. 

The  stream  thrust  function  can  then  be  simplified  for  axial  thrust  in  the  follow¬ 
ing  equation: 

7 ZT 

Sat  =  Ui(l  +  (3.4) 

The  value  of  Sa*,  where  i  represents  a  given  index  along  the  x-axis,  was  then  plotted 
for  all  values  of  x.  The  difference  between  the  exit  and  inlet  conditions  served  as  a 
measure  of  total  stream  thrust  produced  in  the  nozzle. 

3.3  Efficiency  Computation  Model 

For  this  study,  the  efficiency  of  MGD  acceleration  was  described  by  the  following 
relationship, 
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(3.5) 


APout 

11  =  aF" 

APout,  is  the  increase  in  power  due  to  the  MGD  accelerator.  This  can  be 
determined  by  subtracting  a  baseline  (no-MGD)  power  output  from  any  investigated 
MGD  case.  The  power  output  due  to  the  MGD  accelerator  can  be  expressed  in  terms 
of  the  increase  in  the  product  of  thrust  and  velocity [17].  The  velocity  for  this  study 
will  be  taken  as  the  change  in  velocity  between  the  inlet  and  exit  of  the  nozzle,  A u. 

F  F 

A Pout  [( X- )  hi  Alt]  MGD  [(  •  )  Baseline  (3-6) 

m  m 

where  rh  is  the  mass  flow  rate  which  is  defined  as  puA  in  which  A  is  the  area,  of  the 
cross  section  of  the  nozzle.  Since  m  remains  constant  by  definition,  it  did  not  matter 
where  it  was  calculated.  For  consistency,  it  was  calculated  at  the  exit  plane. 

Assuming  that  the  magnetic  field  originates  from  a  permanent  magnet,  which  is 
not  augmented  with  any  electromagnetic  fields,  power  input,  A Pin,  is  actually  the  sum 
of  the  power  requirements  of  the  two  major  components  that  make  MGD  acceleration 
possible.  These  included  the  power  requirement  that  maintained  the  electric  potential, 
A PinA4,,  and  that  which  maintained  the  conductivity  pattern,  A Pin<T. 

A^Pin  =  ^  APinA0ii  +  A  Pina  (3.7) 

where  n  is  the  index  for  the  electrode.  In  this  study,  there  were  only  three  electrodes 
so  n  ranged  from  1  to  3. 

The  value  of  A Pina  could  be  found  using  a  separate  solver  for  ionization  based 
on  flow  parameters[21].  However,  This  study  does  not  attempt  to  calculate  the  power 
requirements  due  to  the  conductivity  pattern.  It  should  just  be  noted  that  the  effi¬ 
ciency  values  obtained  are  inflated  due  to  the  omission  of  this  term.  In  subsequent 
chapters,  this  inflated  efficiency  will  be  referred  to  as  pseudo  efficiency. 
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The  power  requirements  for  the  electric  potential,  A 0,  whose  SI  units  are 
Volts[6],  were  determined  from  basic  principles  in  electricity [30].  Figure  3.1  is  a  dia¬ 
gram  of  the  electrodes.  Each  set  of  electrodes  consisted  of  an  anode  and  a  cathode 
based  on  electric  potential,  0.  The  space  between  was  a  medium  of  air  with  its  own 
conductivity  value,  a.  The  distance  between  the  electrodes  was  l. 


Figure  3.1  Diagram  of  electrodes. 


For  simplicity,  E  was  assumed  to  be  constant  in  the  z-direction, 


A(j)  =  (j)b  -  4>a  = 


f*b  ^ 

E-ds  =  E  dx  =  El 

i  Jo 


(3.8) 


Given  the  above  equation,  Eqn.  3.8,  and  Ohm’s  law,  as  defined  in  the  rest 
frame  from  Eqn.  2.27,  Ohm’s  law  can  be  rewritten  in  terms  of  the  width  of  the  nozzle 
between  the  electrodes,  £,  and  the  electric  potential,  A0  in  Eqn.  3.9. 


3  =  »■[#  -  UB] 


(3.9) 
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Here,  the  — UB  term  results  from  the  property  that  U  x  B  is  working  against  the 
electric  potential.  For  reference,  charge  flux,  j,  is  defined  as, 


3 


1 

A 


(3.10) 


where  I  is  equal  to  current  and  A  is  the  cross  sectional  area  of  the  electrode.  Assuming 
E  is  in  the  z-direction,  the  terms  in  Eqn.  3.9  are  rearranged  to  solve  for  /  in  Eqn.  3.11. 


I 


UB] 


Since  electric  power  is  defined  as, 


(3.11) 


P  =  A  <t>I  (3.12) 

Eqn.  3.11  and  Eqn.  3.12  can  be  used  to  determine  the  power  input  requirements  due 
to  the  electric  potential. 


A  Pin*,  =  oA[AEl  _  u  B(A<p)}  (3.13) 

Eqn.  3.13  was  then  used  to  determine  the  value  of  PinArj,  in  Eqn.  3.7. 

3-4  Computational  Domain 

The  computational  domain  for  the  MGD  accelerator  was  derived  from  an  AJAX 
type  MGD  scramjet  engine  investigated  by  Gaitonde  [10].  The  dimensions  of  the 
MGD  accelerator  portion  of  the  nozzle  appear  in  Figure  3.2.  The  measurements  were 
provided  by  Harrington [16].  The  inlet  area,  is  0.1075  m2  and  the  exit  area,  is  0.4335 
m2.  This  difference  was  on  account  of  a  4°  divergence  in  the  top  and  each  side  wall. 
The  exit  to  inlet  area,  ratio  was  4.032. 
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Figure  3.2  Three  dimensional  MGD  accelerator  geometry. 

The  grid  used  for  the  computations  was  built  from  Harrington’s  [16]  domain  in 
which  Gridgen  vl5  was  used  to  develop  a  laminar  grid.  This  laminar  grid  has  ten 
cells  in  the  boundary  layer.  The  boundary  layer  thickness,  determined  using  laminar 
boundary  layer  theory,  was  5.38  mm.  The  orthogonal  grid  solver  was  used  to  mesh 
the  outer  grid  lines  perpendicular  to  the  walls.  The  index  dimensions  of  the  grid  were 
120  x  46  x  56  nodes  which  provided  297,  000  cclls[16]. 

3.5  Computational  Parameters 

The  flow  conditions  used  at  the  accelerator  inlet  were  set  using  the  Mach  8 
design  flight  conditions  which  were  employed  by  Harrington[10].  Table  3.1  lists  the 
free  stream  conditions  and  the  inlet  conditions  used  in  the  numerical  computations. 
For  reference,  Table  3.1  also  lists  the  non-dimensional  values  of  the  inlet  conditions 
when  ratio  of  specific  heat  is  set  to  1.23.  [16] 

Throughout  the  results,  dimensional  values  are  computed  from  non-dimensional  data 
using  the  equations  listed  in  Table  3.2.  These  values  are  the  fixed  parameters  in  the 
scientific  method. 
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Table  3.1  Freestream  and  MGD  Accelerator  Inlet  Conditions. 


Free-stream 

Conditions 

Inlet 

Conditions 

Non-Dimensional 
Inlet  Conditions 

Velocity 

2380.45  m/s 

1938.2  m/s 

0.8142 

Temperature 

250  K 

2892.0  K 

11.568 

Pressure 

1285.76  Pa 

143116  Pa 

1.4101 

Mach 

8 

1.9148 

1.9148 

Table  3.2  Dimensional  Conversion  Factors. 


Multiply  by 

Multiply  by 

Symbol 

Expression 

Value 

Units 

Position 

X 

u 

0.6 

m 

Y 

Z 

Velocity 

U 

U0 

2380.45 

m 

s 

V 

w 

Pressure 

p 

Po 

1285.76 

Pa 

Density 

p 

Pa  _  „ 

287.0-Tn  1  0 

0.01792 

kq 

Temperature 

T 

T0 

250 

K 

Magnetic  Flux  Density 

By 

Bo 

10 

T 

Electric  Potential 

0 

UqLoBo  =  </>o 

14282.7 

Volts 

Conductivity 

a 

PO-Uo  _  T 
BlL0 

0.710961 

mho/m 
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This  study  investigates  the  effect  of  two  independent  parameters.  The  first  of 
which  is  the  conductivity  pattern.  There  were  two  assumed  conductivity  patterns 
which  were  designed  to  simulate  electron  beam  ionization.  In  both  patterns,  the  non- 
dimensional  a*  is  set  to  unity  in  the  center  of  the  accelerator  and  it  decreases  in  a 
modified  Gaussian  distribution  to  a  given  value  at  the  walls[13].  Conductivity  pattern 
2  has  a  larger  area  of  high  conductivity  in  the  axial  direction. 

Figure  3.3  and  Figure  3.4  show  the  conductivity  patters  on  two  cut  planes  in 
the  accelerator  nozzle  to  illustrate  the  three  dimensional  nature  of  the  conductivity. 
Conductivity  pattern  1  shown  in  Figure  3.3  has  a  smaller  range  in  the  x-direction 
than  conductivity  pattern  2  in  Figure  3.4.  Both  conductivity  patterns  assume  that 
(To  =  0.710961  A  Z  —  0  m  slice,  and  a  Y  —  0.15  m  slice  were  chosen  to  convey 
the  pattern  along  the  axial  coordinates  of  accelerator. 


Y 


■  0.7 
0.6 
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□  0.1 
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Figure  3.3  Three  dimensional  conductivity  pattern  1. 


The  dimensional  values  of  a  are  used  to  determine  the  magnetic  flux  density, 
Bq,  which  may  vary  nominally  from  2  Tesla  to  10  Tesla[16].  By  setting  the  interaction 
parameter,  Q,  to  unity  and  rearranging  the  terms  Eqn.  2.47  becomes: 
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Figure  3.4  Three  dimensional  conductivity  pattern  2. 


B0  = 


(3.14) 


in  which  Lo  is  0.6,  and  po  and  Vo  are  the  free-stream  values  from  Table  3.2.  This  means 
that  the  values  of  Oo  and  Bo  depend  on  one  another  such  that  if  a  cro  of  0.710961 
is  chosen  then  the  accelerator  requires  a  magnetic  flux  density  of  10  T;  or  if  a  <Jo  of 
17.778  is  chosen  then  2  T[16].  For  our  study,  B0  was  assumed  to  be  10  T. 

The  second  independent  parameter  in  this  study  was  the  load  factor,  /C,  which 
was  given  in  Eqn.  4.1.  The  magnetic  flux  density  is  non-dimensionalized  to  unity. 
This  means  that  changes  in  the  load  factor  cause  changes  in  E0,  the  magnitude  of  the 
applied  electric  field.  Thus  the  load  factor  becomes  a  measure  of  the  applied  electric 
field  intensity  [16]. 

The  Lorentz  force  given  as  j  x  B  is  applied  to  accelerate  the  flow  in  the  MGD 
solution  runs.  The  magnetic  flux  density,  B ,  is  oriented  in  the  positive  y  direction. 
For  the  load  factors  greater  than  1.0,  an  electric  current  is  produced  in  the  negative 


3-9 


z  direction  which  results  in  an  accelerating  Lorentz  force  in  the  positive  x  direction. 
For  the  load  factors  less  than  1.0,  the  induced  effects  are  greater  than  the  electric  field 
so  electric  energy  is  generated  causing  decelerated  flow  in  the  x  direction. 

3. 6  Computational  Cases 

The  flow  solution  was  allowed  to  converge  without  any  MGD  interaction.  This 
case  was  called  the  Baseline,  and  its  converged  solution  was  used  as  the  initial  condi¬ 
tions  for  all  of  the  subsequent  cases  [16]. 

The  subsequent  cases  involved  MGD  interaction.  MGD  acceleration  is  defined 
by  a  load  factor,  /C,  greater  than  1.0  and  generation  is  defined  by  a  load  factor  less 
than  1.0.  The  independent  parameters  of  conductivity  pattern  was  also  varied. 

The  generating  load  factors  investigated  were  0.8  and  0.9  for  both  sigma  pat¬ 
terns.  The  accelerating  load  factors  ranged  from  1.2  to  1.6.  Table  3.3  shows  the 
naming  convention  used  in  Chapters  4  and  5  for  the  investigated  cases. 

Table  3.3  Naming  Scheme  for  Investigated  Cases. 

Load  Factor  Conductivity  Conductivity 


(/C) 

cr-pattern  1 

cr-pattern  2 

0.8 

K08 

K08_SIG2 

0.9 

K09 

K09_SIG2 

1.2 

K12 

K12J3IG2 

1.4 

K14 

K14.SIG2 

1.6 

K16 

K16.SIG2 

The  K10  case  still  has  a  voltage  potential  across  the  electrodes  except  the  electric 
held,  E,  is  set  equal  to  the  UB  term.  K10  is  thus  different  from  the  baseline  case 
which  does  not  include  electrodes  at  all. 

The  voltage  potential  for  each  electrode,  A 0  was  determined  using  the  following 
equation: 


A  0  =  YiiUrefB 


(3.15) 
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where  £  is  the  width  between  the  two  electrodes. 
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4.  Results 


The  MGD  accelerator  described  in  Chapter  3  was  used  in  this  numerical  study  to 
investigate  the  influence  of  accelerating  and  decelerating  load  factors.  Table  3.3  shows 
the  nomenclature  used  for  the  different  runs.  The  Non-MGD  baseline  case  was  run  to 
convergence  and  then  used  as  a  reference  in  later  MGD  solutions.  The  objective  was  to 
determine  performance  outcome  based  on  the  specific  thrust  and  the  efficiency  of  the 
nozzle.  The  first  section  provides  the  specific  thrust  and  quasi  1-D  interpretation  of 
the  stream  thrust  function.  This  serves  as  an  overview  for  the  following  three  sections 
which  interpret  the  behavior  of  the  3-D  flow  fields.  Then  non-MGD  baseline,  the 
K16  MGD  acceleration,  and  K08  MGD  deceleration  cases  are  investigated  in  detail  to 
show  characteristics  of  the  flow  fields.  The  detail  flow  field  figures  of  all  other  MGD 
cases  investigated  are  presented  in  Appendix  A.  Last,  the  efficiencies  of  the  MGD 
accelerated  cases  are  presented. 

4-1  Specific  Thrust 

Specific  thrust  for  each  case  was  calculated  using  Eqn.  3.1.  The  results  were 
arranged  in  the  bar  graph  shown  in  Figure  4.1.  Cases  K12,  K14,  and  K16  increased 
thrust  as  expected  for  both  conductivity  patterns.  These  cases  accelerated  the  flow 
which  caused  an  increase  specific  thrust  over  the  baseline  case.  The  generators  decel¬ 
erated  the  flow  which  caused  a  decrease  in  thrust. 

In  order  to  interpret  the  specific  thrust  results,  a  quasi  1-D  plot  of  stream 
thrust,  Sa,  was  investigated  over  the  nozzle  length  in  Figure  4.2  and  Figure  4.3  for 
conductivity  patterns  1  and  2.  In  both  of  these  figures,  the  K12,  K14,  and  K16 
cases  showed  increase  in  thrust  respectively  across  each  electrode,  interspaced  with 
moderate  increases  in  thrust  due  to  the  expansion  nozzle.  In  general,  the  specific 
thrust  in  SIG2  was  greater  than  that  of  SIG1  because  the  greater  conductivity  allowed 
for  better  current  flow. 

The  behavior  of  the  generator  stream  thrusts  at  the  first  electrode  is  of  inter¬ 
est.  Although,  the  generator  cases  ultimately  decelerated  the  flow  as  expected,  their 
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Figure  4.1  Specific  Thrusts 
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Figure  4.2  Stream  Thrust  Function  for  Conductivity  1 
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Figure  4.3  Stream  Thrust  Function  for  Conductivity  2 


stream  functions  increased  at  the  first  electrode.  The  explanation  for  this  behavior  is 
simple.  The  load  factors  were  set  using  the  reference  velocities  of  the  baseline  case. 
When  A0  in  Eqn.  3.15  is  fixed,  the  load  factor  changes  with  the  local  velocity,  Uref, 
as  shown  in  this  rearrangement  of  the  equation  for  voltage  potential: 


A  0 

iUrefB 


(4.1) 


The  field  of  electric  potential  is  not  limited  to  the  region  directly  between  the  elec¬ 
trodes.  In  fact,  the  leaking  field  causes  the  flow  to  decelerate  before  reaching  the 
electrode.  When  the  flow  has  reached  the  electrode,  it  has  already  decelerated  slightly 
from  the  baseline  case.  The  deceleration  continues  in  a  gradual  manner,  due  to  the 
weak  ionization  at  the  beginning  of  the  electrode.  This  sets  up  the  conditions  for 
a  change  in  load  factor  behavior  as  the  flow  decelerates  beyond  the  point  at  which 
UB  is  no  longer  greater  than  E  (see  Eqn.  4.1).  The  spike  in  the  stream  function  is 
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characteristic  of  the  Gaussian  distribution  for  a,  since  the  electrode  becomes  more 
effective  as  the  load  factor  decreases. 

The  rest  of  the  generator  cases  follow  a  path  of  acceleration  followed  by  imme¬ 
diate  deceleration  between  the  first  and  second  electrodes.  This  occurred  because  the 
potential  across  the  second  set  of  electrodes  was  set  using  the  local  velocity  of  the 
baseline  case.  Since  acceleration  occurred  across  the  first  electrode,  the  local  velocity 
was  greater  than  expected  and  ended  up  artificially  decreasing  the  load  factor  to  that 
of  an  MGD  generator. 

In  general,  SIG2  had  greater  specific  thrusts  for  the  accelerator  and  less  spe¬ 
cific  thrusts  for  the  generators  when  compared  to  the  results  for  SIG1.  The  percent 
increases  in  specific  thrust  are  shown  in  Table  4.1  below. 


Table  4.1  Percent  Difference  in  Specific  Thrust,  from  baseline  case. 

Load  Factor  Conductivity  Conductivity 
/C  Pattern  1  Pattern  2 


0.8 

-31% 

-23% 

0.9 

-22% 

-26% 

1.2 

27% 

30% 

1.4 

54% 

60% 

1.6 

84% 

95% 

As  expected,  thrust  decreases  for  the  /C  =  0.8  and  /C  =  0.9  cases  for  both 
conductivity  patterns.  The  magnitude  of  the  percent  change  is  greater  for  the  /C  =  0.8 
case  than  the  /C  =  0.9  case  in  the  SIG1  pattern.  This  is  expected  because  there  should 
be  more  deceleration  for  the  smaller  load  factors.  However,  the  opposite  is  true  for  the 
SIG2  pattern.  It  is  uncertain  what  caused  this  increase  in  percent  difference;  so  more 
work  needs  to  be  done  in  the  area  of  the  effect  of  conductivity  on  MGD  generators. 
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4-2  Baseline  Flow  field 


The  baseline  solution  for  the  conditions  and  geometry  described  in  Chapter  3 
was  computed  without  MGD  interaction.  The  Mach  number  in  the  x-direction  of  the 
baseline  case  is  shown  in  Figure  4.4.  Although  no  load  factor  is  set  in  the  baseline  case, 
the  cut  planes  at  x  —  0.72647, 1.67647,2.62647  are  at  the  midlines  of  the  electrodes. 
The  flow  accelerates  from  Mach  1.96  at  the  exit  of  the  combustor  to  about  Mach  3.06 
(2480  m/s)  at  the  exit  of  the  nozzle  as  expected. 
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Figure  4.4  Mach  Number  Contours  for  Baseline  Case 


Figure  4.5  shows  the  baseline  temperature  contours.  The  flow  temperature 
decreases  as  the  flow  expands,  as  expected.  The  temperature  of  the  fluid  is  2,  658  K 
in  the  middle  of  the  inlet  plane;  and  1,763  K  at  the  exit  plane.  Higher  temperature 
regions,  3,  677  K,  are  located  on  the  wall  surfaces  due  to  stagnation. 

Baseline  pressure  contours  are  shown  in  Figure  4.6.  The  pressure  decreases  as 
expected  due  to  geometry  expansion  in  the  nozzle.  The  pressure  at  the  inlet  reaches 
a  maximum  of  162,537  Pa  throughout  the  entire  inlet  plane  including  the  walls.  The 
average  exit  pressure  is  21,  711  Pa  in  the  exit  plane. 
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Figure  4.5  Temperature  Contours  for  Baseline  Case 
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Figure  4.6  Pressure  Contours  for  Baseline  Case 
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4-3  Acceleration  Case,  /C  =  1.6 

Flow  field  acceleration  occurs  when  the  electric  and  magnetic  fields  are  ori¬ 
ented  such  that  the  resulting  Lorentz  force  accelerates  the  flow.  The  K12,  K14,  K16, 
SIG27K12,  SIG2_K14,  and  SIG2_K16  acceleration  flow  fields  were  very  similar.  There¬ 
fore,  only  the  K16  case  is  presented  in  the  results.  The  flow  held  plots  for  the  rest  of 
the  cases  can  be  found  in  Appendix  A. 

The  Lorentz  force  vectors  for  K16  case  on  the  midplane  in  the  y  and  z  directions 
are  in  Figure  4.7.  The  forces  act  in  an  accelerating  direction  throughout  the  entire 
electrode  regions  as  seen  in  the  x-z  midplane  on  top.  The  vectors  accelerate  the 
how  inward  toward  the  center  at  the  leading  edge  of  the  electrodes  and  change  their 
direction  along  the  electrodes  to  be  parallel  with  the  how.  At  the  end  of  the  electrodes, 
the  vectors  push  the  how  toward  the  wall.  The  x-y  midplane  force  vectors  show  that 
the  largest  accelerating  forces  are  at  the  locations  of  maximum  conductivity,  near  the 
center  of  the  nozzle  in  the  x  direction. 


Figure  4.7  K16  Lorentz  Force  Vectors  (top:  x-z  midplane,  bottom:  x-y  midplane) 
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The  Lorentz  force  has  the  effect  of  increasing  the  flow  velocity  and  mach  number. 
Figure  4.8  shows  how  the  mach  number  increases  across  each  electrode.  The  largest 
increase  in  Mach  number  occurs  in  the  center  of  the  flow  field  with  respect  to  the 
x-z  plane.  This  is  very  similar  to  the  baseline  case.  However,  the  accelerated  case 
features  long  regions  of  slower  Mach  number  along  the  walls.  This  is  probably  due 
to  a  combination  of  the  Lorentz  forces  pushing  against  the  wall  and  increases  in  local 
temperature  along  the  wall.  The  maximum  Mach  number  at  the  exit  of  the  nozzle  is 
3.13,  which  is  a  2.5%  increase  over  baseline  exit  Mach  number.  Negating  the  effect 
temperature  has  on  Mach  number,  the  velocity  at  the  exit  is  2,730  m/s  which  is  a 
10%  increase  in  the  exit  velocity.  The  Lorentz  force  has  accelerated  the  flow. 


Figure  4.8  K16  Mach  Number  Contours 

The  electric  field,  created  by  the  load  factor  /C  =  1.6,  is  shown  in  Figure  4.9. 
This  vector  field  shows  how  some  of  the  field  lines  in  the  x-z  midplane  travel  from  one 
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set  of  electrodes  to  the  next  set  of  electrodes  on  the  same  wall.  This  is  known  as  Hall 
effect  which  describes  the  phenomenon  that  charged  particles,  like  electrons  and  ions, 
move  in  a  perpendicular  direction  to  an  electric  and  magnetic  field.  [30]  Also,  leakage 
occurs  on  the  inlet  and  exit  of  the  nozzle.  This  is  expected  since  there  is  no  way  to 
insulate  the  center  of  the  flow  field.  The  end  of  the  x-y  midplane  shows  how  some  of 
the  vectors  turn  back  into  the  nozzle. 


_i 
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Figure  4.9  K16  Electric  Field  Vectors  (top:  x-z  midplane,  bottom:  x-y  midplane) 

Figure  4.10  shows  contours  of  the  electric  potential  throughout  the  flow.  It  is 
interesting  to  note  that  the  potential  between  the  second  and  third  set  of  electrodes 
along  the  walls  in  non-zero  as  if  there  were  one  big  cathode  and  anode  pair.  It  seems 
as  if  the  electric  field  created  by  the  electrodes  may  travel  from  the  anode  and  cathode 
in  the  z-direction  as  well  as  to  the  adjacent  electrode  in  the  x-direction.  The  forward 
and  reverse  E  field  vectors  in  Figure  4.9  support  this  explanation.  This  would  also 
account  for  the  Lorentz  forces  in  the  positive  and  negative  z-directions  that  result  fore 
and  aft  the  electrodes  in  Figure  4.7.  However,  this  result  is  more  likely  a  combination 
of  both  the  Hall  effect  and  the  bi-directional  flow  of  the  E  field. 
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Figure  4.10  K16  Electric  Potential  Contours 
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A  side  effect  of  MGD  acceleration  is  the  increase  in  the  temperature  along 
the  walls  due  to  Joule  heating.  Figure  4.11  shows  the  temperature  contours  for  the 
/C  =  1.6.  As  expected,  the  temperature  decreases  in  a  diverging  supersonic  duct. 
However,  the  temperature  increases  by  more  than  an  order  of  magnitude  along  the 
walls  with  the  electrodes.  Figure  4.12  zooms  in  on  the  surface  of  the  wall  at  the  trailing 
edge  of  an  electrode.  Since  the  velocity  is  very  slow  and  the  electric  potential  is  at 
its  maximum  at  the  electrodes,  Joule  heating  causes  abnormally  large  temperatures. 
In  addition,  these  extremely  high  temperatures  are  probably  also  due  to  the  type  of 
boundary  conditions  enforced  along  the  walls.  More  in  depth  analysis  is  needed  on 
the  proper  type  of  boundary  condition  to  employ. 
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Figure  4.11  K16  Temperature  (top:  x-z  midplane,  bottom:  x-y  midplane) 
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Figure  4.12 


Close  up  of  K16  Temperature 
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4-4  Generation,  /C  =  0.8 

Flow  field  generation,  or  deceleration,  occurs  when  the  electric  and  magnetic 
fields  are  oriented  such  that  the  resulting  Lorentz  force  decelerates  the  flow  by  pulling 
energy  out  of  the  flow  field  in  the  form  of  electric  potential.  Since  the  potential  is 
kept  constant  on  the  electrodes,  in  a  real  life  scenario,  the  resulting  energy  would  be 
extracted  as  current.  The  K08,  K09,  K08_SIG2,  and  K09_SIG2  cases  provided  similar 
flow  fields.  As  discussed  in  the  previous  section  on  specific  thrust,  the  generation 
cases  demonstrated  some  tendencies  toward  acceleration  at  the  first  electrode  and 
reverted  to  the  expected  deceleration  behavior  at  the  following  electrode  due  to  the 
conductivity  of  the  flow.  The  K08  case  is  discussed  in  this  section.  The  rest  of  the 
MGD  generation  flow  field  plots  can  be  found  in  Appendix  A. 

The  Lorentz  force  field  vectors  for  K08  case  are  shown  in  Figure  4.13.  The 
largest  deceleration  force  occurs  at  the  second  electrode.  This  happens  because  the 
load  factor  at  the  second  electrode  is  artificially  reduced  by  the  acceleration  at  the  first 
electrode  and  the  conductivity  is  a  maximum  at  the  second  electrode  so  the  Lorentz 
force  is  very  strong  there.  The  third  electrode  shows  a  more  moderate  decelerating 
force  on  account  of  less  conductivity  at  this  electrode. 

The  Mach  number  contours  are  shown  in  Figure  4.14.  The  Mach  number  de¬ 
creases  between  the  second  and  third  electrodes,  and  then  increases  before  it  reaches 
the  exit  plane.  The  second  and  third  electrodes  act  like  MGD  generators,  followed  by 
a  region  of  acceleration  due  to  the  nozzle  geometry  expansion.  This  can  be  confirmed 
by  the  absence  of  vectors  behind  the  third  electrode  shown  in  Figure  4.13.  This  is 
also  supported  by  the  plot  of  specific  thrusts  shown  in  Figure  4.2. 

Figure  4.15  shows  the  electrical  field  vectors.  The  electrical  field  around  the 
first  electrode  looks  similar  to  that  of  the  K16  acceleration  case  as  compared  to  Fig¬ 
ure  4.2.  There  is  a  noticeable  change  in  the  vector  field  around  the  second  and  third 
electrodes:  characterized  by  a  reverse  in  the  flow  field  direction  which  causes  a  decel¬ 
erating  Lorentz  force. 


4-13 


Figure  4.13  K08  Lorentz  Force  Vectors  (top:  x-z  midplane,  x-y  midplane) 


Figure  4.14  K08  Mach  Number  Contours 
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X 

Figure  4.15  K08  Electrical  Field  Vectors  (top:  x-z  midplane,  x-y  midplane) 

The  pressure  contours  are  shown  in  Figure  4.16.  The  pocket  of  relatively  lower 
pressure  at  the  second  set  of  electrodes  is  most  likely  caused  by  the  interface  between 
the  acceleration  and  expected  MGD  generation.  Over  the  later  portion  of  the  nozzle, 
the  decrease  in  pressure  is  retarded,  as  expected,  by  the  decelerating  body  forces. 

The  temperature  contours  are  shown  in  Figure  4.17.  It  is  interesting  to  note 
that  the  temperature  is  not  as  hot  at  the  walls  as  in  the  K16  case.  However,  the 
temperature  of  the  flow  held  in  the  middle  of  the  nozzle  is  much  hotter  than  that  of 
the  K16  case.  This  is  due  to  the  slowing  down  of  the  how. 
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Figure  4.16  K08  Pressure  Contours 
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Figure  4.17  K08  Temperature  Contours  (top:  x-z  midplane,  bottom:  x-y  midplane) 
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4-5  Efficiency 


The  efficiency  results  were  determined  using  Eqn.  3.5  and  Eqn.  3.7  by  setting 
Pina  equal  to  zero.  This  method  lead  to  skewed  results  on  account  of  the  electron 
beam  power  requirement  not  being  included.  These  results  are  referred  to  as  pseudo 
efficiencies  since  they  do  not  include  all  of  the  power  requirements.  Figure  4.18  shows 
the  pseudo  efficiencies. 


Figure  4.18  Accelerator  Pseudo  Efficiency 


Given  Eqn.  3.6  a  greater  specific  thrust,  equates  to  a  larger  increase  in 
power  above  that  of  the  baseline  case.  As  a  result  when  electron  beam  requirements 
are  neglected,  it  would  be  expected  for  the  case  with  the  largest  specific  thrust  to 
also  have  the  largest  efficiency.  However,  this  is  not  the  case  in  Figure  4.18  in  which 
the  K16  case  has  the  smallest  pseudo  efficiency.  This  is  partly  due  to  the  fact  that 
the  voltage  power  requirements  PinA0  was  dwarfed  by  Pout  for  the  lower  load  factors. 
It  is  also  assumed  that  the  power  requirements  for  the  conductivity  pattern  are  the 
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dominant  power  requirements.  This  means  that  even  an  estimation  of  the  efficiency 
would  require  the  Pina  term  in  order  to  make  sense. 

In  order  to  account  for  the  effect  that  electron  beam  power  requirements  may 
have  on  the  efficiency,  a  projected  efficiency  was  calculated.  Figure  4.19  shows  the  ef¬ 
ficiency  that  would  result  given  an  arbitrary  PiUcr .  As  expected,  the  projected  efficien¬ 
cies  show  that  the  higher  specific  thrust  producing  load  factors  have  better  efficiencies 
than  the  lower  ones. 


Figure  4.19  Accelerator  Projected  Efficiency 


It  is  interesting  to  note  that  an  optimistic  efficiency  of  97%  results  if  the  electron 
beam  power  requirement  for  the  second  conductivity  pattern  does  not  exceed  0.2 
gigawatt  when  /C  =  1.6.  However,  more  research  would  have  to  be  done  on  the  power 
requirements  for  an  electron  beam  capable  of  producing  the  conductivity  pattern 
in  this  study.  This  leads  to  the  notion  that  MGD  acceleration  can  be  a  practical 
application,  as  long  as  the  practical  barriers  such  as  material  limitations  are  overcome. 
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5.  Conclusions 


5. 1  This  Study 

The  specific  thrust,  flow  conditions,  and  efficiency  results  for  the  scramjet  nozzle 
accelerator  were  examined  computationally  for  a  laminar  flow  with  a  free  stream 
Mach  number  of  8.  The  changing  parameters  included  load  factors  and  conductivity 
patterns.  The  ultimate  goal  was  to  determine  which  parameters  yielded  the  best 
performance  in  terms  of  specific  thrust  and  efficiency. 

The  efficiency  of  the  MGD  acceleration  could  be  reasonable  if  the  electron  beam 
ionization  was  added  to  the  calculation.  The  pseudo  efficiency  results,  ranging  from  a 
factor  of  5.04  to  7.69,  from  this  study  only  accounted  for  power  requirements  due  to  the 
electric  potential  requirements  across  the  electrodes.  Since  the  power  requirements  for 
electron  beam  ionization  plays  a  significant  role  in  efficiency,  the  pseudo  efficiencies 
alone  where  not  enough  to  draw  conclusions.  Figure  4.19  was  created  to  project 
efficiencies  that  would  result  given  possible  electron  beam  power  requirements.  Given 
a  beam  ionization  power  requirement  of  1  gigawatt  or  less,  the  efficiency  of  the  /C  =  1.6 
case  for  SIG2  could  yield  decent  results.  However,  more  research  would  be  needed  to 
determine  power  requirements  for  each  flow  field. 

A  more  definite  performance  measure  used  in  this  study  was  specific  thrust.  As 
shown  in  Table  4.1,  specific  thrust  increased  for  the  MGD  acceleration  cases,  and 
decreased  for  the  MGD  deceleration  cases.  The  second  conductivity  pattern  with 
more  ionization  yielded  more  acceleration.  The  best  acceleration  was  found  for  the 
/C  =  1.6  case  with  the  second  conductivity  pattern  which  had  a  95%  increase  in 
specific  thrust  over  the  baseline. 

More  research  should  be  done  into  deceleration  load  factors  in  the  nozzle.  The 
/C  =  0.9  case  for  conductivity  2  yielded  a  —26%  difference  in  specific  thrust  which 
fell  outside  of  the  trend  considering  that  the  /C  =  0.8  only  yielded  a  —23%  difference. 
This  may  be  clue  to  the  flow  field  decelerating  enough  to  drift  the  load  factor  in 
Eqn.  ??,  which  has  the  result  of  accelerating  the  flow.  This  effect  may  not  have  been 
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as  noticeable  in  the  SIG1  case  simply  because  the  conductivity  decreased  enough  to 
have  a  only  minute  effect  on  the  flow  field  when  the  load  factor  increased. 

On  a  similar  note,  the  initial  acceleration  at  the  first  electrode  for  all  the  genera¬ 
tion  cases  a  was  due  to  the  fact  that  the  converged  baseline  case  was  used  to  determine 
the  load  factor,  /C,  at  each  electrode.  The  first  electrode  acted  like  an  accelerator  for 
all  the  cases.  Then  the  flow  shifted  to  generation  as  intended.  Eqn.  4.1  accounts 
for  this  effect.  The  same  principle  has  a  similar  effect  on  the  acceleration  cases.  As 
the  flow  accelerated,  the  load  factors  at  each  successive  electrode  became  less  than 
intended.  This  had  the  effect  of  weaker  increases  in  the  stream  thrust  function,  and 
ultimately  less  specific  thrust. 

In  order  to  prevent  the  load  factor  drift,  a  controller  should  be  implemented  to 
vary  the  voltage  with  the  reference  velocity  to  regulate  the  change  in  load  factor.  In 
computational  studies,  this  only  means  recalculating  Eqn.  4.1  based  on  local  param¬ 
eters  throughout  each  successive  iteration.  However,  in  real  world  applications,  this 
would  mean  designing  a  control  system  with  sensors  capable  of  determining  velocity 
of  fluids  without  effecting  the  flow.  Possible  solutions  could  include  lasers,  or  even 
feedback  data  from  a  power  supply.  Since  the  Lorentz  force  is  a  product  of  j  x  B ,  and 
the  load  factor  represents  energy  interaction  with  the  electrodes,  controlling  the  load 
factor  could  be  achieved  through  the  use  of  a  current  pumping  power  supply.  Such 
power  supplies  are  used  in  a  variety  of  applications  such  as  arc  welding. 

An  increase  in  load  factor  had  the  effect  of  increasing  both  specific  thrust  and 
efficiency  for  each  of  the  two  conductivity  patterns.  The  only  drawback  with  increas¬ 
ing  load  factor  is  the  addition  of  Joule  heating  at  the  electrodes.  SIG2,  the  larger 
conductivity  pattern,  also  yielded  higher  specific  thrusts.  However,  efficiency  con¬ 
clusions  could  not  be  drawn  concerning  the  change  in  conductivity  pattern  since  the 
power  requirements  for  conductivity  were  undetermined.  MGD  control  may  have  a 
plausible  future. 
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5. 2  Future  Work 


The  accuracy  of  the  computational  code  when  determining  Joule  heating  at  the 
walls  remains  uncertain.  A  study  could  be  designed  to  test  the  limitations  of  the 
computational  code.  It  might  be  interesting  to  find  if  there  are  certain  load  factors 
that  cause  problems. 

More  investigation  could  also  be  done  into  the  efficiency  of  the  generating  con¬ 
ditions.  Macheret  et  al.  of  Princeton  university  developed  a  method  of  determining 
efficiency  of  the  generator  using  computational  parameters [23].  He  also  designed  a 
computational  study  that  may  have  been  more  efficient  by  placing  small  scale,  on 
the  order  of  millimeters,  electron  beam  windows  directly  in  front  of  each  of  ten  elec¬ 
trodes.  It  might  be  interesting  to  compare  the  ionization  results  of  this  study  to  those 
of  Macheret  et  al.  Their  study  had  efficiencies  on  the  order  of  13%  with  electron  beam 
power  requirements  of  8.72  kilowatts [23].  Since  this  study  is  on  a  completely  different 
scale,  comparison  of  those  results  is  not  possible.  However,  it  would  be  interesting  to 
set  up  a  computational  experiment  that  could  compare  the  two  methods. 
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Appendix  A.  Other  Cases 


Figure  A.l  K09  j  x  B 
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Figure  A. 2  K09  Mach  Number 
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Figure  A. 3  K12  j  x  B 


Figure  A. 4  K12  Mach  Number 
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Figure  A. 5  K14  j  x  B 


Figure  A. 6  K14  Mach  Number 
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Figure  A. 8  K08_SIG2  Mach  Number 
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Figure  A. 9  K09J3IG2  j  x  B 


Figure  A.  10  K09_SIG2  Mach  Number 
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Figure  A.  11  K12J3IG2  j  x  B 


Figure  A.  12  K12_SIG2  Mach  Number 
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Figure  A.  13  K14_SIG2  j  x  B 


Figure  A.  14  K14_SIG2  Mach  Number 
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Figure  A.  15  K1623IG2  j  x  B 
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Figure  A.  16  K16_SIG2  Mach  Number 
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Figure  A. 17  SIG1  Velocity 
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Figure  A.  18  SIG2  Velocity 
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Figure  A.  19 
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Figure  A. 20  SIG2  Temperature 
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Figure  A. 21  SIG1  Pressure 
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Figure  A. 22  SIG2  Pressure 
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